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ABSTRACT 

We construct the first "face-on" Tully-Fisher (TF) relation for 24 galaxies 
with inclinations between 16° and 41°. The enabling measurements are integral- 
field, echelle spectroscopy from the WIYN 3.5m telescope, which yield accurate 
kinematic estimates of disk inclination to ~ 15°. Kinematic inclinations are 
of sufficient accuracy that our measured TF scatter of 0.42 mag is comparable 
to other surveys even without internal-absorption corrections. Three of four 
galaxies with significant kinematic and photometric asymmetries also have the 
largest deviations from our TF relation, suggesting that asymmetries make an 
important contribution to TF scatter. By measuring inclinations below 40°, we 
establish a direct path to linking this scatter to the unprojected structure of disks 
and making non-degenerate dynamical mass-decompositions of spiral galaxies. 

Subject headings: galaxies: fundamental parameters — galaxies: structure - 
galaxies: kinematics and dynamics — galaxies: spiral 



1. Introduction 

One common characteristic of Tully-Fisher (TF; Tully & Fisher 1977) studies of the 
relation between luminosity and rotation speed is the selection of spiral galaxies with in- 
clinations (z) above 45°, and typically with i ~ 60°. The inability to determine accurate 
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and precise inclinations for galaxies with i < 40° has limited studies attempting to link 
structural parameters, such as lopsidedness and ellipticity, to TF scatter (Franx & de Zeeuw 
1992, Zaritsky & Rix 1997). Likewise, combining vertical stellar velocity dispersions and disk 
rotation speeds has been thwarted for lack of precise inclinations for nearly face-on galaxies 
(e.g., Bottema 1997). Hence, non-degenerate dynamical mass-decompositions of galaxies 
have not yet been made. 

One solution to the above problems is to use nearly face-on galaxies to study the TF 
relation. Here the detailed planar structure of the disk can be studied unprojected, while 
the vertical component of the disk stellar velocity dispersion is favorably projected. The 
problem is now observational: How can accurate and precise inclinations be measured for 
small inclination angles? 

For most TF studies, inclinations are calculated from photometric axis-ratios, assum- 
ing circular disks, and correcting for a constant, intrinsic flattening (e.g. Tully & Fouque 
1985; Courteau 1996). Bars, spiral structure, and intrinsic ellipticity cause systematic er- 
rors in photometric inclination and position angles (PA), dominate the apparent ellipticity 
at low inclination, and skew apparent inclinations to larger values (Andersen et al. 2001; 
Courteau et al. 2003). Hence there is a triple penalty at low inclinations for TF studies: 
(1) Erroneously large inclinations lead to under- correcting the velocity projection. (2) Ran- 
dom errors diverge since photometric disk axis ratio measurements have constant errors as 
a function of axis ratio. (3) For long-slit studies, PA-mismatch will always yield systematic 
underestimates of the projected rotation speed. 

These pitfalls can be mitigated if kinematic inclinations are derived from two-dimensional 
velocity fields (e.g., Schommer et al. 1993; Verheijen 2001) because disk scale heights are 
irrelevant, and velocity fields between 1 and 3 scale lengths are typically less aberrated than 
the light distribution by spiral structure, bars or warps (Courteau et al. 2003; Briggs 1990). 
A limit of ~ 40° had been set on the minimum derived kinematic inclinations for nearby 
galaxies from HI maps (Begeman 1989), but this minimum value is not fundamental; it is 
notionally adopted as such. 

We have discovered that kinematic inclinations, derived from high-quality, optical kine- 
matic data collected with integral field units, are sufficiently accurate and precise to construct 
a TF relation at inclinations as low as 15°. In this letter we present an analysis of the first 
TF relation for a sample of nearly face-on galaxies. 
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2. The Photometric and Kinematic Data Set 

We obtained R and J-band images and Ha velocity fields of 39 spiral galaxies meeting 
the criteria of Andersen et al. (2001): Galaxies appear nearly face-on (i < 30°), are of 
intermediate type (Sab-Sd), with no clear signs of interactions, bars or rings, and are in 
regions of low Galactic extinction. Our sample and observations are described in Andersen 
(2001). 

CCD images were acquired over 17 nights between May, 1998 and January, 2001 at 
the WIYN 3.5m 2 , the KPNO 2.1m, and the McDonald Observatory 2.7m telescopes. Total 
magnitudes are defined as asymptotic fluxes, as measured from growth curves, calibrated 
on the Kron-Cousins system using Landolt (1983) standards observed on four photometric 
nights. Zeropoints for other nights are internally boot-strapped, where possible, from repeat 
observations. We use Hq = 70 km s _1 Mpc^ 1 and correct magnitudes for Galactic absorption 
(Schlegel, Finkbeiner & Davis 1998). (Virgo in-fall affects recession velocities by <1%; 
Paturel et al. 1997.) We use the mean R— I color of 0.52±0.08 mag for our calibrated sample 
to estimate the R magnitude for several galaxies with only accurate J-band photometry, and 
add the 0.08 mag range in R — I in quadrature as an additional uncertainty for these sources. 

Ha velocity fields were obtained at an instrumental resolution of -£r ~ 13, 000 using 
the DensePak integral field unit (Barden, Sawyer & Honeycutt 1998) on the WIYN 3.5m 
telescope over 11 nights between May, 1998 and January, 2000. [PGC 56010, a low surface- 
brightness galaxy, was re-observed using SparsePak (Bershady et al. 2003) in May, 2001.] 
The Ha emission lines were resolved, with typical signal-to-noise levels of ~ 35 and cen- 
troid precision of ~ 1.5 km s _1 . DensePak observations yield velocity centroids for ~ 140 
independent positions within each galaxy out to 2.2 to 5.8 disk scale lengths (typically 3.3). 

3. Measuring Inclinations Below 40° 

The most critical measurement required to construct a "face-on" TF (FOTF) relation 
is the derivation of accurate and precise disk inclinations. As we have argued, there are 
fundamental reasons why two-dimensional kinematic maps provide the only viable data for 
measuring these inclinations. Using a method similar to that described in Andersen et al. 
(2001), we fit the DensePak data with a single, hyperbolic-tangent velocity-field model. 
From these fits, we derive measurements of kinematic inclinations, rotation speeds (V vot ), 
and recession velocities for our sample. 
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Of the 39 galaxies in our sample we measure reliable deprojected rotation velocities, 
V circ = V rot / sini, for 26 (see Table 1; two galaxies lacking calibrated R and /-band mag- 
nitudes are excluded). This subsample has V rot < 120 km s _1 with inclinations and errors 
from 16°1^ to 41°ljj. Ten additional galaxies have measured, deprojected rotation veloci- 
ties, but kinematic inclinations consistent with zero within 68% confidence intervals. Their 
mean inclination is 8°. Three remaining galaxies have no projected circular velocity field; 
they appear truly face-on. Our results paint a consistent picture of being able to measure 
inclinations reliably down to i ~ 15°. 

Why can we measure accurate inclinations for nearly face-on galaxies? For disks with 
inclinations between 20° and 30°, the peak difference in projected velocity is only 2 to 3 
km s _1 between disks with inclinations differing by 5°. To measure this difference requires 
< 2 km s _1 centroiding precision per fiber, or many samples of the velocity field at a precision 
at least comparable to this difference. The latter, which we achieve, is preferable because 
of random motions. We find typically 3 km s~ x rms variations between our model and the 
data (Andersen 2001). 

Historically, HI data have been claimed to be of insufficient quality for deriving galaxy 
inclinations below 40° (Begeman 1989). One aspect of the problem has been that tilted 
ring models include a weighting scheme which, when typically implemented, minimizes the 
importance of velocity measurements at azimuthal angles which are most influenced by 
changes in inclination. For example, Begeman's simulations use a cos# weighting scheme, 
where 9 is the angle from the major axis in the galaxy plane. This weighting is done to 
minimize the effects of warps, but in so doing over half of the signal is thrown out. For the 
same reason, data are often not considered beyond 9 of 30° or 45°, throwing out yet more of 
the signal. 

Our DensePak observations sample the inner disks of galaxies, where warps tend to be 
negligible. Hence we can use the entire velocity field to maximize the influence of inclination 
in the model velocity fields. Figure 1 shows deprojected position-velocity diagrams for two 
representative galaxies in different bins of 9. Excellent agreement between azimuthal bins 
indicates a good model fit and lack of systematics. We find bidimensional echelle spectroscopy 
to be of sufficient quality to yield deprojected rotation velocities, V TOt / sini at precisions 
better than 10% at inclinations above 28°; a 20% precision on V rot / sini is reached for 
inclinations of ~ 20°. 
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4. FOTF vs Conventional TF Relations 

To evaluate our FOTF data, we compare to the TF samples of Courteau (1997) and 
Verheijen (2001). The salient difference between our sample versus Courteau's and Verhei- 
jen's is that their galaxies have % > 40° and % ~ 65°. Courteau's sample shares comparable 
scale lengths, colors, surface brightnesses, distances and uses spatially-resolved Ha long-slit 
emission as the kinematic tracer. We use the cleaner "Quiet Hubble Flow" (QHF) subset of 
Courteau's sample and convert Courteau's Lick r-band photometry to Cousins .R-band using 
it! — r — —0.36 (Fukugita et al. 1995). Verheijen collected multi-band (including Cousins 
i?-band) photometry and HI synthesis maps for a volume-limited sample from the nearby 
Ursa Major cluster. We select the Ursa Major subsample with M R < —19, corresponding to 
the faint end of the FOTF and Courteau samples, and adopt the HST Key Project distance 
(Sakai et al. 2000). Both Courteau and Verheijen use an internal-absorption correction of 
the form A { = 7 log a/b which minimizes the apparent TF scatter. (We have subtracted from 
Courteau's data his zeropoint correction of 0.40 mag, which corresponds to his absorption 
correction at i — 70°.) Since our sample is nearly face-on, we apply no internal- absorption 
corrections to our photometry. 

An inspection of the three samples in Mr versus log(2U rot / sin 2) of Figure 2 highlights 
the differences and similarities in their TF relations. There is excellent agreement between 
the Courteau and Verheijen TF slopes and zeropoints. There is also an excellent agreement in 
the zeropoint of the FOTF with Verheijen and Courteau TF relations near L* at M R ~ —22, 
but the shallower slope of the FOTF appears striking. However, slopes vary between —7 
and —9 for both the Courteau and Verheijen samples depending on the method of velocity 
measurement, the extinction correction prescription and the sample selection. Verheijen 
(2001) notes the above absorption corrections steepen TF slopes. However, these corrections 
do not alter the slope of the FOTF relation. While empirical corrections (e.g., Giovanelli 
et al. 1994, Tully et al. 1998) work well for galaxies with 40° < i < 80°, the corrections 
for galaxies with i < 30° are small (A R ~ 0.1 mag). Since in practice one can systematize 
velocity measurements, a larger FOTF spanning a broader range in luminosity may provide 
a new test of the veracity of TF-based internal-extinction corrections. 

When we restrict our sample by luminosity or V TOt , we find slopes between -5 and -8; 
a FOTF subsample with 2.3 < log(2V r circ ) < 2.7 shares the same slopes and zeropoints as 
the Courteau and Verheijen culled samples (Figure 1). This implies that only the extreme 
rotators, fast and slow, differ. This may be due to luminosity-dependent extinction correc- 
tions discussed above. It is more difficult to present a plausible systematic error whereby the 
high luminosity objects rotate faster and the low luminosity objects rotate slower, and the 
intermediate objects are unaffected. However, with only six galaxies in these regimes, the 
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discrepancy between slopes may be influenced by our small-number statistics coupled with 
relatively large errors on log V/ sini 

If the FOTF is to be a useful tool for studying internal extinction or, e.g., correlations 
of galaxy properties to TF scatter, the impact of inclination measurements on TF residuals 
needs to be examined. Figure 3 shows the relation between TF residuals and inclination 
along with the predicted envelope of TF residuals introduced by inclination error. We find 
the scatter for our sample is comparable to the scatter in the culled Courteau and Verheijen 
samples: 0.42 mag for ours versus 0.38 mag for the Verheijen and Courteau samples used 
here, each about their respective regressions (Figure 2). Six galaxies labeled in Figures 2 
and 3 are noted in Table 1 as either being photometrically or kinematically asymmetric, or 
having small recession velocities (hence uncertain distances). Inclination errors alone aren't 
the cause of their large TF residuals. Removing these six galaxies, the scatter for our face-on 
sample about its shallow regression drops significantly to 0.3 mag, and is consistent with our 
measurement errors. 



5. SUMMARY AND DISCUSSION 

We have presented a TF relationship for 24 galaxies with a mean inclination of 26° - the 
first time the TF relation has been measured for a sample of galaxies with inclinations less 
than 40°. Accurate and precise kinematic inclinations are the key to studying TF relation at 
low inclination, and velocity fields are the essential ingredient for their measurement. We are 
able to construct this FOTF relation from integral-field Ha velocity fields for four primary 
reasons: (1) DensePak and SparsePak are efficient for gathering Ha velocity fields with 
moderately high spatial resolution and medium spectral resolution at high signal to noise. 
(2) Ha velocity fields generally extend out to ~ 3 disk scale-lengths so that warping, which 
dominates the outer radii of some galaxy disks, is unimportant. (3) Instead of using multiple 
tilted rings, we successfully fit galaxy disks with a single velocity-field model, simultaneously 
using data at all radii and phase angles. (4) We equally-weight all data, making our fits 
sensitive to the signal constraining inclination, which peaks in the velocity field at ~ 45° 
from the major axis. 

The FOTF relation measured from our data has a zeropoint and scatter comparable 
to conventional TF relations but possibly a shallower slope. The bulk of the FOTF scatter 
is consistent with measurement errors, as found by Verheijen (2001). Of the four galaxies 
which show large kinematic and photometric asymmetries, three fall well off the TF relation. 
While our statistics are small, this is consistent with the claim by Barton et al. (2000) that 
kinematic asymmetries are an important source of TF scatter. 
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The FOTF can play an important role in understanding the nature and interpretation of 
the TF slope, zeropoint and scatter in three ways. First, internal extinction is less important 
for face-on galaxies minimizing inclination-dependent effects. For inclined TF samples, if 
correlated errors exist between de-reddened luminosities and de-projected rotation-speeds, 
TF scatter will be increased. If extinction depends on luminosity or velocity (mass), the 
corrected TF slopes may differ from their true values. 

Second, measurements of asymmetry, ellipticity, and lopsidedness ideally are made at 
low inclinations. Several studies have suggested non-axisymmetric dark halo mass distribu- 
tions such as lopsidedness (Zaritsky & Rix 1997; Swaters et al. 1999) or ellipticity (Franx 
& de Zeeuw 1992; Andersen et al. 2001) could be sources of TF scatter. However, without 
access to inclination measurements for face-on systems, these claims could not be directly 
evaluated. Our small FOTF sample suggests such links are possible, and show that such 
studies can now be pursued wholesale. 

Finally, the mass-to-light ratio of spiral disks is an essential ingredient for interpreting 
the TF zeropoint, slope, and scatter. Only limits can be placed on mass decompositions (Bell 
& de Jong 2001) without direct dynamical estimates of total and disk mass. This can be 
accomplished via rotation curves and vertical stellar velocity dispersions of disks. The latter 
are best measured in face-on systems, and hence the FOTF unlocks the door for determining 
the fundamental mass-budget of spiral galaxies. 

We thank E. Bell, S. Courteau, and B. Tully for helpful comments. We acknowledge 
NSF AST-9970780 and AST-0307417. 
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Table I. Tully-Fisher Parameters 





% 


^circ 


Mr 


PGC 


(deg) 


(km s _1 ) 


(mag) 



2162 
3512 
5673 
6855 a 
7826 b 
14564 a 
16274 
23333 
23598 
24788 a 
26140 
28310 
31159 a 
32638 
33465 
36925 
38908 
39728 6 
46767 
56010 
57931 
58410 
72144 
72453 



25.811* 

30.1127 

21.01^ 

38.810 
40.710 

+2.9 
-3.1 
+2.3 
-2.7 
+3.8 
-4.5 



30.5 
30.6 
26.9 



16 2 +5 ' 2 



28.2l|° 

qcr q+2.8 

+3.0 
-4.0 
+3.7 
-4.9 

18.61^ 



22.6 
22.0 



22.6; 



r3.5 

-4.5 



31.8; 

22.7- 



26.01^ 

h3.3 
-3.6 
h2.6 
-3.0 

26.7l« 3 3 6 
21.51^ 
29.01^ 
27.91^ 

ni O+4.0 



oi o+38.5 

164.31™ 
160.6™ 

84.3l£° 

49.0l^ 3 2 
181.81^ 
222.41^ 
133.01™ 

ic9 o+189.8 
iuz '°-100.7 



41.51^ 112.81 2 . 7 ; 4 



qi 7 q+88.1 

oi i -o_ 61 g 

110.91}™ 
220. 1™ 

234.25l^;4 
309.8l^J 
195.2l 5 7 
211.01^ 

- L ^- z, -26.8 
9qq 9+43.4 

ycr O+70.7 
< d -<J-43.2 

192.41^8 
231.4l|;? 
110 6+ 42 - 8 

1±U.U_ 30 g 

212.4l^ 5 3 6 



-20.42±0.03 
-21.17±0.03 
-20.89±0.05 
-20.64±0.05 
-19.30±0.06 
-21.48±0.05 
-22.03±0.05 
-20.67±0.06 
-21.69±0.05 
-21.70±0.05 
-22.86±0.05 
-20.38±0.07 
-21.19±0.07 
-21.80±0.05 
-22.41±0.05 
-22.06±0.07 
-21.96±0.06 
-20.32±0.06 
-22.79±0.07 
-19.97±0.07 
-21.83±0.05 
-22.11±0.10 
-21.27±0.03 
-22.23±0.03 



a Significant photometric and kinematic 
asymmetry. 

^ Heliocentric < 3000 km S -1 . 
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Fig. 1. — Azimuthally-deprojected position-velocity diagram for PGC 23598 (Top) and PGC 
6855 (Bottom). Each datum is a separate fiber measurement. Velocities, color coded as 
receding (red) or approaching (blue), are marked as open circles for 9 < 20° and filled 
squares for 20° < 9 < 60°. Solid curves trace the projected model rotation curve along the 
major axis, x 2 rnaps of the 68% and 90% confidence limits on model inclination and rotation 
velocity are inset. PGC 6855 is photometrically lopsided and kinematically disturbed; it falls 
well off the TF relation (Figure 2). PGC 23598 has the lowest sample inclination but falls 
directly on the TF relation. 
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Fig. 2. — Tully-Fisher relations for our sample of galaxies with i = 26° (FOTF, white and 
black) and more-inclined samples of Courteau (1997; red) and Verheijen (2001; blue). FOTF 
"outliers" (filled black circles) are defined in text and Table 1. The heavy solid line with 
slope -4.90 is the best fit to the FOTF data (with or without outliers) using a backward 
regression which allows for intrinsic scatter and heteroscedastic errors (BCES[X|Y]; Akritas 
& Bershady 1996). The heavy dashed line is the best fit to the FOTF sample between 
2.3 < logVJ-ot < 2.7; it is indistinguishable from the best fits to Courteau and Verheijen 
(light solid line with slope -8.2). 
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Fig. 3. — TF residuals versus inclination for sample galaxies. 68% confidence limits on each 
galaxy reflect uncertainty in kinematic inclination and rotation velocity propagated through 
the best-fit TF relation. The envelope defined by the two, solid curves enclose the area in 
which 68% of points are expected to fall if all TF scatter were due solely to uncertainties in 
kinematic inclinations. "Outlier" galaxies (see text and Table 1) are labeled. That 5 of the 
6 labeled galaxies fall outside the envelope suggests that factors besides inclination errors 
make these galaxies TF outliers. 



